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' Abstract 

Nuclear effects in J/tjj hadro- and electroproduction on nuclei are considered in framework 
of reggeon approach. It is shown that screening regime which holds for electroproduction at 
x f ^ 0.7 and for hadroproduction at xf > —(0.3 -r- 0.4) is changed with anti-screening regime 
^ ■ for smaller xf values. 

m . 

^r) . Heavy-quark production on nuclei provides an important information on strong-interaction mech- 

anisms and is intensively discussed (see e.g. refs ^ |21 E] for review of experimental data, discussion 
of some phenomenological models and more references) . In this note we will discuss a phenomenon of 
£f) , changing screening regime to anti-screening one in J/ ift hadro- and electroproduction when x^ = xp 

decreases in the framework of the BCKT model jlj |Sj based on the reggeon approach. For J/tp 
• hadroproduction it happens in the region of negative . An investigation of this region can provide 

an important information on dynamics of charmonium production and may discriminate different 
Qh! dynamical models [2J |S] . 

Nuclear effects are usually discussed in terms of conventional power-law parameterization Fa{x^) oc 
Fn(xip) ■ A a ( x ^\ where Fn(x^) ( Fa(x^) ) is inclusive J /if) production cross section on a nucleon 
(on a nucleus). The function a(x^) characterizes nuclear effects at different longitudinal momentum 
^ | fraction x^. 

Experimental data on J/ip hadroproduction reveal a striking contradiction with simplest theo- 
retical expectations. Experimentally jZHHUH] function a{x^) decreases from 0.93-^0.95 at ~ to 
values ~ 0.75 at x^ ~ 0.8 thus indicating an increase of the absorption as x$ increases. Formation- 
time mechanism predicts an opposite behavior. Instead of expected scaling with pfj^ experiment 
shows an approximate Feynman scaling with x^ [B] . Comparison of J/ ip production data at different 
energies exhibits an explicit breakdown of QCD factorization theorem for this process |10j . 

On the other hand when comparing x dependence of a for J/tp hadroproduction with one for 
light-quark particle hadroproduction (inclusive production of pions, nucleons, lambdas, etc.), one 
observes the same trend - a small absorption at x ~ 0, large absorption at x — 1 and approximate 
Feynman scaling. The difference is only quantitative one - a is about 0.4 -i- 0.5 at x ~ 1 for light 
quarks instead of ~ 0.7 for heavy quarks. This behavior of a(x) for light hadrons allows a natural 
explanation in Regge theory. Small absorption at x ~ is due to Abramovskii-Gribov-Kancheli 
(AGK) cancellation for inclusive spectra and increasing absorption at high x is due to violation 
of AGK rules because of momentum conservation requirement j!2). 
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In ref.jl] the model for heavy quark and lepton pair production was constructed taking into 
account the most essential aspects of reggeon approach. In the spirit of parton picture it 
was suggested that a fast projectile looks due to quantum-mechanical fluctuations like a cloud 
of virtual particles consisting of light partons (quarks and gluons, or light-quark hadrons, mainly 
pions) and (with some small probability) of heavy partons (say cc pair). 1 Different constituents of 
this fluctuation interact with nuclear matter and determine a dependence upon atomic number A. 
There are both elastic and inelastic interactions which from the viewpoint of reggeon diagrams are 
the different discontinuities of the same reggeon diagram. As a result there exist definite numerical 
relations between the discontinuities of different types (AGK rules jllj). 

As it was stressed in ref.jl] one has to distinguish two types of reggeon diagram cuttings. For 
diagrams of the first class the registered particle is contained inside a pomeron and appears in 
intermediate state only due to cutting of this pomeron. In this case the AGK rules are always valid 
and all contributions of many-pomeron diagrams to inclusive spectra cancel. It happens because 
each additional pomeron can be both cut and uncut, and due to the opposite signs of inelastic and 
elastic rescatterings they cancel. As a result, only one-pomeron diagram contribution survives giving 
a linear A dependence of spectra. This situation is typical for low x particles. 

Another situation occurs if a particle to be registered is contained inside "vertex" of reggeon 
diagram (we use a term "vertex" for part of diagram common for several pomerons attached to it). 
In this case AGK cancellation is not valid in general. The striking manifestation of this mechanism 
comes from A dependence of inclusive spectra at large x close to 1 (momentum conservation mecha- 
nism \1'2\). The fastest particles belong to the "vertex" but not to the pomeron, and corresponding 
discontinuities give a ~ 1/3 for large A compared to a ~ 1 for small x. 

Similar effect exists if one observes the "vertex" particle of some particular type, e.g. the J/ip 
meson. There are two sources of nuclear effects in this case J3]. The first one is connected with 
rescatterings of light partons in fluctuation on different target nucleons or, alternatively in the anti- 
laboratory frame, with fusion of fluctuations originated from different nucleons of fast nucleus. This 
mechanism plays a small role at present energies pi a b < 10 3 GeV [I] and its contribution is inessential 
for anti-screening effects we are interested in. 

Another contribution is connected to rescatterings of the charmed state itself. One of the impor- 
tant ingredient of the model of ref. ,3] is an assumption that it is not the J/ tp meson that propagates 
along nucleus but some primary colorless system containing both cc quarks (in color state) and light 
quarks (to screen this color charge). This state denoted as X can be DD mesons, or D*D* state 
or something else, but the crucial point is that this system is of large size and therefore strongly 
interacts with nucleons (with cross section of order of 20 -r- 30 mb). Only at last stage this system 
X converts into J/ijj meson (or in tp') which is registered experimentally. 

After inelastic rescattering the state X can be lost for registration (e.g. due to conversion to 
a state with small projection into J/ip). One can describe this effect quantitatively introducing 
a probability e to find the state X after its inelastic rescattering. Then for e = 1 one has exact 
AGK cancellation, for e = the situation is similar to momentum conservation mechanism for light 
mesons at x close to 1 (no AGK cancellation) with a ~ 1/3, and for < e < 1 the situation is 
an intermediate one with incomplete AGK cancellation. Moreover, the momentum of the state X 
decreases after each inelastic rescattering. This redistribution of the longitudinal momentum of the 

1 In reggeon theory a mechanism of emitting these partons is supposed to be a multiperipheral one. Creation of 
heavy quarks in this model can be considered as an example of "intrinsic charm" mechanism which was discussed by 
Brodsky et al. |14II1U) . In reggeon approach one can consider on equal footing any heavy component (e.g. "an intrinsic 
lepton pair" and so on). 
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state X is essential for anti-screening effects discussed below. 

Let us denote the momentum fraction distribution of system X immediately after its production 
as Fq(x). The state X has some projection into state J/ij) (as well as into ip', etc.) 2 . The x^ 
distribution can be obtained from Fq(x) by convolution with some projection function G^: 



where the notation was used 



or, in the rapidity variables, 



F { J ) \x f ) = F 1 ®G^ (1) 



1 dz 

!/\ ,/)(,•)- / T f(z)g(x/z) , (2) 



(/®</)(v) = ^ <%f(t)g(y-0 , (3) 
J o 

where y = Y — y is a rapidity in the anti-lab frame. 

The state X travelling through nuclear matter suffers both elastic and inelastic rescatterings. 

As in ref. [I] the change of x distribution due to inelastic rescattering will be described as a convo- 

l 

lution with some function eG(z), where G(z) is normalized (J G{z)dz/z = 1), and the parameter e 

o 

determines a probability to have again the state X after inelastic interaction (though with different 
momentum). After k inelastic rescatterings of the state X we have fe-fold convolution: 

F fc W = Fx ® eG® - • -®eG ® G^ , (4) 

s v y 

k times 

It is essential that the operation of convolution is commutative and associative, so we can first 
convolute F\ with G^ in the Eq. @ and thus get the fc-fold convolution of the J / tp spectrum on the 
nucleon, Ffo(x), with rescattering functions G: 

Ff = F$® eG®- - -®eG . (5) 

k times 

In order to obtain the J/ip inclusive spectrum on a nucleus, the functions F^ should be weighted 
with cross sections a A for k inelastic rescattering after X production: 



oo 

Ft = Y. F ^)°A +l) ■ (6) 

fc=0 



The explicit form of the weights cr^ +1 ^ depends upon energy region |3]. However the change 
is rather smooth and for simplicity we will use the same low-energy expressions as in ref. [I]: 

r r T A(b) 

°A +1) = <N J d2b J Q dv eM-vax)(va x ) k /kl , (7) 



2 For simplicity we suggest that only one state X has a noticeable projection into J/ip. In order to describe properly 
the observed difference of nuclear dependencies for J/ip and tp' , one has to consider at least two primary states, Xi 
and Xi, with different projections to J/ip and ip' and slightly different interactions with nuclear matter. 
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where aZl is the cross section of J/ t/j production on a nucleon, ax is the cross section of interaction 
of the state X with a nucleon, and Ta(6) is two-dimensional nuclear density as a function of impact 
parameter b (nuclear profile). These formulas correspond to rescatterings ordered in longitudinal 
direction (low-energy regime) and after summation over k give well known optical-type formula for 
the production cross section. Thus we come to Eq.© for J/tj) nuclear inclusive cross section where 
Ffc(x) and are defined by Eqs.© and 0. 

Let us discuss parameterizations of distribution functions entering the model. The function 
F^ N {x^) can be determined from experimental data on J/ip production on nucleon. Its form strongly 
depends on a type of a projectile a. For hadron beam it can be parameterized as 

!&(*) = C7.(l-*)* , (8) 

where n ~ 2 for pion and (3 p ~ 3 for proton beam. For J/i/j production by photons (and electrons) 
the form of this function is strongly different: produced ifi particles are concentrated near ~ 1, 
and Xih distribution can be parameterized in a power-like way 

f{}W=|5/, (9) 

with /3 7 > 6. We choose the value /3 7 = 6. 

As to the function G(x), it is natural to suggest its similarity to the function F^(x). For 
simplicity we choose exactly the same form 

G(x) = fix , (10) 

with (3 considered as a free parameter. 

Terms of the series © in number of rescatterings k decrease fastly with k, and it is enough to take 
into account several first terms (4-j-5). The corresponding convolutions at fixed impact parameter b 
are estimated analytically, and the b integration was performed numerically. We use for calculations 
the Woods-Saxon parameterization of nuclear density with standard values of parameters taken from 
ref.[El. 

Before comparison of theoretical predictions with experimental data let us discuss firstly quan- 
titative character of absorptive corrections. As we discussed above at small x ~ we have AGK 
cancellation, and the extent of its violation is proportional to 1 — e, a probability to lose the state 
X after inelastic rescattering. The particle can be lost for registration not only if it disappears 
after interaction but also if it loses part of its momentum and is shifted into different region in x. 
As a result at large x values the AGK rules are violated and absorption effects can be observed 
(momentum conservation effect). However it is only redistribution (not a loss) of particles for the 
whole x region (see Fig.l for illustration of this fact). For cross section integrated over x screening 
(absorption) effects are connected only with 1 — e value but not with momentum conservation. For 
J/tj) photoproduction the absorption effects were usually analyzed only for integrated cross section 
and this is a reason for rather small values of effective cross sections extracted for J/ifi absorption 
(3 -r- 5 mb). We want to demonstrate here 3 that as a function of x nuclear effects can be rather 
big (with effective cross section of order of 20 -r 30 mb) but integrally they cancelled because of 
cancelling contributions from screening and anti-screening regions. 

This is demonstrated at Fig. 2 where we show description of the NMC data ^H] on ratio R(Sn/Be) 
of J/ip electroproduction cross sections as a function of x = xp. Since in this case the J/ip spectrum 

3 First this phenomenon was discussed in the talk based on unpublished results by K. Boreskov and A. Kaidalov. 
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is concentrated at x close to 1, the anti-screening regime takes place already at x < 0.8. Note that 
theoretical predictions at small x value (x < 0.6) are less reliable due to their sensitivity to small- x 
parameterization of functions F^j and G(x). 

For J/ip hadroproduction the anti-screening regime can be expected (see Fig. lb) only for negative 
x. Fig. 3 and Fig. 6 show that for J ftp and ip' hadroproduction the change of regime happens at 
x ~ —(0.3 -T- 0.5). This value depends on parameters of the model - the less are momentum losses 
(the larger is fi) and the less is e the more negative x's are necessary for the anti-screening onset 
(see Fig. 3). For example, for the value of (3 = 4 the crossover happens already at x ~ —0.2. 

It is important to note that it is not entirely adequate to analyze such delicate effects as variation 
of dependence on atomic number with change of x or px in terms of the standard parameterizations 
A a . In general, a is not only a function of x and pt but depends also on atomic number A. To 
illustrate this fact we show at Fig. 4 the ratios of spectra R{A/N) versus A at fixed x^ which by no 
means are described with a power dependence. The same effect is demonstrated in Fig. 3: there is 
a noticeable difference of solid and dashed theoretical curves which correspond to different choices 
of A ranges. Therefore the a value extracted from experiments depends significantly on a chosen 
range of A and on a way of analysis. To exclude this uncertainty it is preferable to compare with 
theory ratios of cross sections for different nuclei. Description of such data from E772 experiment is 
shown in fig. 5. 

A difference in nuclear effects in J/tp and tp' hadroproduction should be connected in the model 
with presence of several primary states, say X and X' , which interact with nuclear matter in different 
ways and have different branchings into J /tp and ip'. In Fig. 6 we give as an example a comparison 
of data on tp' hadroproduction from refs |H] with model calculation with a single X state at the 
same parameter values as in fig. 2. The two-channel generalization is evident. 

Let us mention that in some theoretical models 0111] used for description of J/tp, ip 1 hadroproduc- 
tion in the negative x region the function a(x) decreases (screening effects increase) as \x\ increases 
in this region. 

Thus, an experimental investigation of nuclear effects for charmonium production in the whole 
region of x can provide a valuable information for testing different dynamical models. In the reggeon 
approach one can distinguish three different regions in rapidity or in Feynman variable for J/ip 
hadroproduction - in central region one has strong cancellation of screening diagrams and a is close 
to one (small absorption cross section); at x > 0.2 due to violation of cancellation rule the absorption 
effects increase and effective cross section is about 20 30 mb; at negative x less than —(0.3 -j- 0.4) 
the anti-screening regime should be observed. Data on A dependence of charmonia production in 
the negative x region will be available soon from HERA-B experiment |17j . 
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Figure 1: Change of the original spectrum (solid curves) after 3 (dashed curves) and 5 (dotted 
curves) rescatterings for: (a) photoproduction and (b) hadroproduction of J/ip. 
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Figure 2: Description of data ^S] on ratio R(Sn/C) for J/tp electroproduction. The curve corre- 
sponds to model calculations at parameter values a = 30 mb, e = 0.85, (3 = 12. 
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Figure 3: Description of data Ej on nuclear dependence of J/ip hadroproduction. Curves corre- 
spond to model calculations for different ways of a extraction: solid curve - from the ratio R(Fe/C) 
and dashed curve - from R(W/C). The dotted line extracted from R(W/C) corresponds to e = 0.95. 
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Figure 4: Ratios R(A/D) vs A at fixed x F = 0.26. 
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Figure 5: Ratios R(A/D) vs xp for different atomic numbers A. 
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Figure 6: Description of data |S] on nuclear dependence of tp' hadroproduction. The curve 
corresponds to model calculations at parameter values a = 30 mb, e = 0.85, j3 = 12. 
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